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The physics of quantum walks on graphs is formulated in Hamiltonian language, both for simple 
quantum walks and for composite walks, where extra discrete degrees of freedom live at each node 
of the graph. It is shown how to map between quantum walk Hamiltonians and Hamiltonians for 
qubit systems and quantum circuits; this is done for both a single- and multi-excitation coding, and 
for more general mappings. Specific examples of spin chains, as well as static and dynamic systems 
of qubits, are mapped to quantum walks, and walks on hyperlattices and hypercubes are mapped to 
C ■ various gate systems. We also show how to map a quantum circuit performing the quantum Fourier 

transform, the key element of Shor's algorithm, to a quantum walk system doing the same. The 
results herein are an essential preliminary to a Hamiltonian formulation of quantum walks in which 
coupling to a dynamic quantum environment is included. 
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In many quantum-mechanical systems at low energies, the Hilbert space truncates to the point where the system 
is moving between a set of discrete states (which may however be very large in number) . In this case we can describe 
the system, with complete generality, as equivalent to a system in which a particle (which may itself possess internal 
degrees of freedom) 'hops' between a set of 'nodes', or 'sites', on some graph - the nodes of this graph can then be 
^ . identified with states in the Hilbert space of the original system. 

The hopping amplitudes between nodes are just the transition amplitudes in the original Hamiltonian, so that the 
topology of the graph is entirely determined by these transition amplitudes. In general we may allow the Hamiltonian 
to be time-dependent, so that both the hopping amplitudes and the on-site energies are allowed to change. We can 
also allow the internal state of the hopping particle to couple to its coordinate on the graph. 

In path integral language, one can think of the trajectory of a quantum particle moving between 2 nodes A and B 
on this graph as a 'quantum walk', made up of a succession of discrete hops. The amplitude to go from A to B is 
^J*. then given by summing over all possible paths (or 'walks') between them, with the appropriate amplitudes. 

Formulated in this way, the problem of a 'quantum walk' is very familiar to most physicists, and has in fact been 
• i-h , under study since the very beginning of quantum mechanics. Notable examples come from solid-state physics (where 
^ 1 particles hop around both crystalline lattices [l[ and disordered systems of various topology Q), from quantum 
magnetism [3[ (where an assembly of spins makes transitions between different discrete spin states), from atomic 
physics and quantum optics (where one deals with discrete atomic states, and where in the last few years 'optical 
lattices' have come under study (4J), and from a large variety of problems on different sorts of graph in quantum 
statistical mechanics 

Quantum Walks and Quantum Information: A certain class of quantum walks has recently come under study in the 
context of quantum information processing Q . These walks are intended to describe the time evolution of quantum 
algorithms, including the Grover search algorithm and Shor's algorithm. The general idea is that each graph node 
represents a state in the system Hilbert space, and the system then walks in 'information space'. In some cases explicit 
mappings have been given between the Hamiltonian of a quantum computer built from spin-1/2 'qubits' and gates, 
and that for a quantum particle moving on some graph [6|, [7| . More generally, the mapping between a walk and an 
algorithm is most transparent for spatial search algorithms with the local structure of the database. 

The quantum dynamics between two sites A and B on a given graph has been shown for certain graphs to be 
much faster (sometimes exponentially faster) than for a classical walk on the same graph [TH, [3] • It has also been 
argued that quantum walks may generate new kinds of quantum algorithm, which have proved very hard to find. 
Those algorithms based on quantum walks proposed so far fall into one of two classes The first is based on 

exponentially faster hitting times @, H, H, G3 > where the hitting time is defined as the mean 'first passage' time taken 
to reach a given target node from some initial state. While several examples have been found, such as the 'glued-trees' 
of Childs et al. [8j| , there is presently no application of these to solve some useful computational problem. The second 
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class uses a quantum walk search [T(| [IH, [l3| providing a quadratic speed-up. In the case of a spatial search, the 
quantum walk algorithms can perform more efficiently than the usual quantum searches based on Grover's algorithm. 
Amongst the graphs so far studied for quantum walks are 'decision trees' @, H, Q and hvpercubesfiol. : quantum 
walks on some other graphs, and their connection to algorithms, were recently reviewed[6|. 

Several recent papers have also proposed experimental implementations of quantum walks for quantum information 
processing [Tfil. in various systems such as ion traps, optical lattices and optical cavities. Some of these proposals 
involve walks in real space, whereas others are purely computational walks (eg., a walk in the Hilbert space of a 
quantum register To our knowledge, two quantum walk experiments have been carried out: a quantum walk 
on the line, using photons [l8j], and a walk on a TV = 4 length cycle, using a 3 qubit NMR quantum computer 
However many experiments over the years, particularly in solid-state physics, have also been implicitly testing features 
of quantum walks. 

The variety of walks that one may consider is quite enormous - one may vary the topology of the graphs, and, 
as we will see below, even quite simple walks may have a complicated Hamiltonian structure on these graphs. Even 
the solid-state and statistical physics literature has only considered a small part of the available graph structures. 
In the quantum information literature, the discussion of walks has so far been confined to a very restricted class of 
graphs and Hamiltonians on these graphs. Attention has focussed almost exclusively on either regular hypercubic 
lattices, on trees (or trees connected by random links), and on 'coin-tossing' walks on lines. Often it is not obvious 
how one might implement these walks in some real experiment - clearly one is not going to be building, for example, 
a d-dimensional hyperlattice! Thus one pressing need, which is addressed in considerable detail in the present paper, 
is to give explicit mappings between the kinds of qubit or gate Hamiltonian that one is interested in practise, and 
quantum walk Hamiltonians. 

Quantum Walks and Quantum Environments: The range of possible quantum walk systems becomes even more 
impressive if one notes that any quantum walker will couple to its environment. In general one needs to understand 
what form the couplings will take, and how they will influence the dynamics of the quantum walk. Typically these 
couplings can be formulated in terms of 'oscillator bath' [2(| HH or 'spin bath' [13, [23[ models of the environment; in 
the case of quantum walks we will see that various couplings to these are allowed by the symmetries of the problem. 
It has been common in the quantum information literature, at least until very recently, to model decoherence sources 
and environmental effects using simple noise sources (usually Markovian) . Results derived from such models are highly 
misleading - they miss all the non-local effects in space and time which result when a set of quantum systems are 
coupled to a real environment, and also give a physically unrealistic description of how decoherence occurs in many 
systems. 

Thus another pressing need is to set up a Hamiltonian description of quantum walkers coupled to the main kinds of 
environment which do exist in Nature, showing how these Hamiltonians transform when one maps between quantum 
walk systems and qubit or quantum gate systems. This then allows a bridge to real experiments. This is actually a 
rather substantial task which is undertaken in a separate paper (24T | . 

Plan of paper. The main goal of the present paper is to set up a Hamiltonian description of quantum walk systems, 
and to give a detailed derivation of the mappings that can be made between quantum walk systems and more standard 
qubit and gate systems. The results are in some cases quite complex, and in order to make them both useful and easier 
to follow we give detailed results for several examples. Two things we do not do in thispaper are (i) incorporate 
couplings to the environment into the discussion - this is the subject of anoth er p aper [24|: and (ii) work out the 
dynamics of walkers for any of the Hamiltonians we derive (see however refs [24| |25| |) . 

In section[H]we begin by setting up a formalism for the discussion of different kinds of quantum walk. In section UTT1 
we then show one may systematically map from different quantum walk Hamiltonians to various qubit systems and 
quantum circuits. This is done first with single- and multi-excitation encoding of walks into many-qubit systems, and 
then more generally; the mappings are illustrated with simple examples, notably walks on a hyperlattice. In section 
IIVI we do the reverse, mapping qubit systems back to quantum walks. This is done first for systems which can be 
maped to spin chains, and then for more general qubit systems, both static and dynamic; to illustrate the mappings 
we discuss various chains and small qubit systems, and show how to map systems implementing the quantum Fourier 
transform to quantum walks. Finally, in the concluding section [V] we summarize our results. 

II. QUANTUM WALK HAMILTONIANS 

In this section we discuss the structure of the different kinds of quantum walk Hamiltonian we will meet. We deal in 
this paper with 'bare' quantum walks (ie., those without any coupling to a background environment). We emphasize 
that in this section (and the next) our primary object of study is the quantum walk, as opposed to, eg., qubit networks 
or quantum circuits. However in section 4 we will be freely mapping between quantum walk systems and other kinds 
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of network. 

We assume, as in the introduction, that the bare walk is defined by the topology of the graph on which the system 
walks, and by the 'on-site' and 'inter-site' terms appearing in the Hamiltonian. We can then begin by distinguishing 
two kinds of bare quantum walk, which we call 'simple' and 'composite', as follows: 

A. Simple Quantum Walk 

The 'simple' quantum walker has no internal states, so that we can describe its dynamics by a Hamiltonian with 
N nodes, each labelled by an integer j € [0, N — 1], of form: 

ij 3 

= -J2^(t)(m\ + \m) + X>(*)Mii (i) 

ij 3 

Here each node j corresponds to the quantum state \ j) — |0), so that \j) denotes the state where the 'particle' is 
located at node j. The two terms correspond to a 'hopping' term with amplitudes Ay(t) between nodes, and on-site 
node energies £j(t), both of which can depend on time. There is no restriction on either the topology of the graph, or 
on the time-dependence of the {Ay (i), ej(t)}. Thus, for example, one can design a pulse sequence for the parameters 
Aij(t) and £j(t), as a method of dynamically controlling the quantum walk. 

Two of the simplest topologies that have been discussed in the literature for quantum walks are d-dimensional hy- 
percubes and hyperlattices. The hypercube simply restricts the simple quantum walk described above to a hypercubic 
graph - its interest resides in the fact that we can map a general Hamiltonian describing a set of d interacting qubits 
to a quantum walk on a d-dimensional hypercube. This mapping is discussed in section ITVl Hyperlattices extend the 
hypercube to an infinite lattice in d dimensions; it is common to assume 'translational symmetry' in the lattice space, 
which means writing a very simple 'band' Hamiltonian 

H=-J2 A o (4^3+^1) = 5>o(p)cpc p (2) 
ij P 

where A Q is a constant, and p is the 'quasi-momentum' (also called the 'crystal momentum' in the solid-state litera- 
ture); the 'band energy' is then 

d 

e (p) = 2A Q ^2 cos(p M a ), (3) 

and the states of the walker can be defined cither in the extended or reduced Brillouin zone of quasi-momentum space. 
In ([3]) we assume a lattice spacing a , the same along each lattice vector; (henceforth we will put a a = 1). All results 
can be scaled appropriately if these restrictions are lifted. 

B. Composite Quantum Walk 

The composite walker has 'internal' degrees of freedom, which can function in various ways. We assume these 
internal modes have a finite Hilbert space, and they can often be used to modify or control the dynamics of the 
walker. Thus we assume a Hamiltonian in which the simple walker couples at each node j to a mode with Hilbert 
space dimension Ij, and on each link {ij} between nodes to a mode with Hilbert space dimension roy, and we have 
a Hamiltonian 

kc = -Y,(Fij(M ij ;t)4d j + H.c) +Y,G j (C i ;t)c\<l j + H Q ({M lj , Cj}). (4) 

ij 3 

This composite Hamiltonian reduces to the simple walker when Fij(Mij]t) — > Ajj(i) and when Gj{Cf,t) —> £j(t). 
We do not at this point specify further what are Fij(Mij\t) and Gj{Cj]t), nor the form of their dynamics (which is 
goverend not only by the coupling to the walker but also by their own intrinsic Hamiltonian H ({M.ij,jCj})), but we 
will study several examples below. The bulk of this paper will be concerned with the simple walker in |T]), which is 
already rather rich in its behaviour. 
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We emphasize that the internal variables are assumed to be part of the system of interest - that is, they are 
not assumed to be part of an 'environment' whose variables are uncontrolled and have to be averaged over in any 
calculation. In the context of quantum information theory these internal variables are assumed to be under the control 
of the operator. For example, Feynman's original model(26j of a quantum computer is a special case of a composite 
quantum walk with Hamiltonian 

H C = -Y.{ F ^-M^+H.c)j : (5) 

ij 

where r corresponds to a set of register spins, where the computation is performed. The walker implements the clock 
of this autonomous computer. Another example of a composite quantum walk is given by the Hamiltonian 

H c = -^^ l S{t-t n )f{Lj;t) (c\c J+ H.c) + H ({Lj}), (6) 

ij n 

in which decisions about where the walker hops to are made at various times U by discrete variables {Lj}. Such 
models include examples where some sequence of pulses acting on the internal walker variables are used to influence 
its dynamics. A simple special case of such Hamiltonians assumes the walk is entirely on a 1-dimensional line, and 
that the discrete variable Lj is just a spin- 1/2 variable - for example, we can assume the form 

= ~\ E E *(* - [(! + *J + (1 - ^cj^cj + H {{fj}), (7) 

j n 

which is just the discrete-time coin tossing Hamiltonian, in which a walker at site j hops to the left/right depending on 
whether the 'coin (ie., spin- 1/2) at this site is up/down, with decisions being made after regular intervals of discrete 
time t Q . Obviously one can cook up many more examples of composite walk systems. 

We have sometimes found it convenient to rewrite both (fTJ) and Q as sums over the original graph Q and an 
ancillary graph Q* formed from the links between the nodes of the graph. Thus we can write, for example, 

Hc = ~ E (MMflt^ej+Hx^+^GjiCirtcfe + H ({M r ,Cj}) (8) 

This representation puts the 'non-diagonal' or 'kinetic' terms on the ancillary lattice on the same footing as the 
'diagonal' or 'potential' terms existing on the original lattice. Such a manouevre can be very useful in studying the 
dynamics of the walker, but we will not need it in this paper. 

In our study in this paper of mappings from quantum walks to systems of qubits and/or quantum gates (or 
vice-versa), we will concentrate on simple walk systems, for two reasons. First, as we will see, the results just for 
simple walks are rather lengthy. Second, a proper discussion of these mappings in a Hamiltonian framework requires 
a treatment of non-local effects in time, which also arise in the discussion of the coupling of the walker to the 
environment. Thus we reserve a detailed treatment of composite walks for another paper. 



III. ENCODING QUANTUM WALKS IN MULTI-QUBIT STATES 

We would now like to map quantum walk systems to a standard quantum computer made from qubits or quantum 
gates. This means that we wish to map from a quantum walk Hamiltonian like |T]), acting on states \j), to a qubit 
Hamiltonian acting on M qubits; and we require an encoding of the node state \j) in terms of the 2 M computational 
basis states. We will use the following notation for the computational basis states, 

\z\z 2 ■ ■ ■ z M ) = \zx) <g> \z 2 ) ® . . . ® \z M ) , (9) 

where Zk S [t, I] (we use spin operators here, instead of the more standard [0, 1], so as to avoid confusion with the 
node indices). 

We now describe two such encodings and the corresponding multi-qubit operators needed to implement the quantum 
walk described by the Hamiltonian J!}, thereby deriving the equivalent qubit Hamiltonian. 
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A. Single-excitation encoding 



Our first encoding implements the quantum walk in an M -dimensional subspace of the full 2 M dimensional Hilbcrt 
space for M qubits. In this sense, this encoding is inefficient in its use of Hilbert space dimension. However, the 
operations can prove to be more easily implementable, requiring only two-qubit terms in the Hamiltonian. 

The subspace we are interested in is spanned by the M-qubit states with only a single excitation - the states with 
only a single qubit in the 'up' state ||) fe state, with all other qubits in the |J.) . state (for all j ^ k). Each node 
of the graph is then encoded via the location of the excitation (in this case, we label the nodes from 1 to N) i.e., 
\k) = ||)x ® |4) 2 O ■ ■ ■ <£> |T)fc ® ■ ■ ■ ® \0)n- l n tms encoding, the general quantum walk Hamiltonian (fl]) is 



2^e,(t)(l+r;), 

3 



(10) 



which consists solely of 2-qubit terms, between each connected pair of qubits, as defined by the graph. This encoding 
allows the implementation of any quantum walk using only two-qubit terms in the Hamiltonian, provided arbitrary 
pairs of qubits can interact. 

To simulate evolution according to Hamiltonian (fTU)) it suffices to be able to explicitly perform controlled evolution 
according to each term in the Hamiltonian. Letting H = Hk, for time-idependent parameters, we use the Trotter 
formula 



-ihtH 



IT 



-ihtH k /N 



N 



(11) 



approaching equality as iV 
unitary 



oo. For time- varying parameters in the Hamiltonian, H{t) 1 evolution is given by the 



U(t, 0) = exp 4 



ih I H{t')d£ 



where exp + is the time-ordered exponential. This can be expanded as the product 

U(t, 0) = U(mS, (to - 1)8) . . . U{8, 0), 



(12) 



(13) 



for small time step 8 = t/m. By choosing 8 sufficiently small, we approximate each term in the Hamiltonian to be 
constant over this time interval, 



U ((n + 1)8, nS) = exp^ 



—iH 



(n+l)S 



nS 



H(t')dt' 



exp [— ihSH (nS)] 



(14) 



Since 8 is small, we then apply the Trotter formula. 

So to simulate the quantum walk on a quantum computer using this single-excitation encoding, we must perform 
unitary operators of the form 

U VJ (e) = e^Mn^+V^ ( 15 ) 
between pairs of qubits representing connected nodes of the corresponding graph, along with the single qubit terms 

t> fe ( e ) = e- ift < (16) 



In this way, this encoding represents a 'physical' walk, of a single spin-up over a network of qubits, defined by the 
pairwise interactions. 

It is interesting to note the scaling of the resources required for such a simulation of a general graph. In terms of 
space, the number of qubits required for a given graph is the corresponding number of nodes. The number of gates 
representing time (assuming only one- and two-qubit operations) is at the very least of the order of the number of 
edges, assuming each qubit is in direct interaction with all others. Details of the scaling of gate resources will depend 
upon the structure of both the graph, and the quantum computing architecture [lj| 
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B. Binary expansion-based encoding 



The most efficient way to encode each node is to use the binary expansion of the integer labelling the node. We start 
from the state at the 'origin' of the quantum walk, and label this state by the ket |0), making this equivalent to the 
qubit 'vacuum state' where all spins are 'down'. Consider, a 2 qubit system. Then we have the mappings |0) = \H), 
|1) = Ut)j 1 2) = ITJ-)) an d |3) = ITT)- The number of qubits required will depend upon the number of nodes of the 
graph - M qubits can encode up to N = 2 nodes. The corresponding many-qubit Hamiltonian for the quantum 
walk depends upon how the nodes of the graph are labelled. We start with the simple example of a free quantum 
walk on the hypercube, before discussing the construction for general graphs, and quantum circuit constructions. 

This encoding represents a walk in information space - the information about the position of the walker is stored in 
a quantum register. A similar construction for the simulation of discrete-time quantum walks on a quantum computer 
was conducted by Fujiwara et al. Results in this section can be viewed as analogous to this work, extended to 

the construction of quantum circuits for simulating continuous-time quantum walks. 

1. Mapping a Hypercube walk to a set of qubits 

Consider hrst the simplest possible quantum walk, where we take ej = (ie., a 'free walk'), and Ay = A D in fl}. 
We also restrict the sum to nearest neighbours, so that H = — A X}<ij> i^l^j + H.c.]. An easily visualised and 
trivial example is a free quantum walk on the regular three dimensional cube. This graph has 8 nodes, so requires 
3 qubits to encode. Figure [T] displays a specific labelling Q and the corresponding qubit encoding To determine the 




FIG. 1: Qubit encoding of a quantum walk on the cubic lattice in three dimensions, using three qubits. 
3-qubit Hamiltonian corresponding to this free quantum walk, one considers a single element, i.e. 

|i)(5| = IUTXUTI 

= liXTI®UXil®lTXT| 



where Pfc = Continuing this process, we obtain 

H = -2A(t + ® I<8>I + I< 
= -4A (t x (g) / (g) i + i c 



ir T , 



H.c 



(17) 
(18) 



which is simply a sum of single qubit terms. 

It is simple to extend this free walk to M-dimensions, where M-qubits are required. Each qubit represents one of the 
M orthogonal directions the quantum walker may move in from each node, and the value of the qubit corresponding 
to that direction gives at which end of that direction the walker is located. The corresponding qubit Hamiltonian for 
the M-dimensional free quantum walk is thus 



(19) 
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The quantum circuit to simulate this Hamiltonian is simply single qubit rotations on each qubit, the angle determined 
by the time of the walk. Scaling of resources for the simulation is trivial - the number of nodes N = log M, while the 
number of gates is the number of qubits, all of which can be applied simultaneously. 

Interactions between qubits are inevitably associated with a 'potential' e,j defined over the nodes, weighted edges, 
and/or next-nearest-neighbour couplings (in section HVl below we derive the relation between the €j and A^- on the 
hypercube and the parameters of a general qubit Hamiltonian) . 



2. General walks and circuit constructions 



From the simple example of the hypercube, we can see how to construct the multi-qubit Hamiltonian corresponding 
to the general quantum walk Hamiltonian using this encoding. Each location/node is now labeled by a bit string 
z = Z\ . . . zm, with"f= 1, J,= 0. A given on-site term in the general quantum walk Hamiltonian |T]) becomes 



M 



M 



A I 



\m=GQ\*k){zk\ = 



(20) 



fe=i 



fc=i 



where T Zk denotes a projection operator. 
For the hopping terms, we have 



c z c w -\- c w c z — 



M M 

\z)(w\ + \w)(z\ = (g) \z k )(w k \ + (g) \w k )(z k \ 

k=l k=l 



(21) 



For each term in the tensor product, either the bit values are equal, and we have a projection operator, or the values 
are opposite, and we have a ladder operator, (t + ,t~), such that 



M 



\z)(w\ + \w)(z\ = ]\(Fl k ) s t'*^5(l-z k -w k )Tjt5(l+z k -w k )Tj-, 

M 

= H(Ft*) s ^-^ (r£ + l (z k - Wk )rlf- S ^- Wk) 



h.c, 



(22) 
(23) 



fe=i 



where 5(x) is the delta function. Expanding the tensor product in terms of Pauli x and y operators, such that the 
addition of the Hermitian conjugate terms ensure only products with even numbers of survive i.e. 



ITITTTIXTITUTI + ITITUTXTITTUI = (jl^t + 



r 4 T 5 T 6 



T 4 T 5 T 6 r 4 T 5 T 6 ) ■ 



(24) 



To simulate the evolution of a general quantum walk on a quantum computer using this encoding, we make use of 
the Trotter formula pip , implying we must be able to implement unitaries corresponding to evolution according to 
each term in the total Hamiltonian. For the onsite/potential terms, this corresponds to unitaries of the form 



U(e) 



_ e -ihe\zlz\ 



(25) 



A simple circuit to implement this unitary |27| uses a single ancilla qubit, initialized in the \[) state, and a multi-qubit 
gate which takes all qubits as input and flips the ancilla qubit if the walker qubits are in the state \z). An example 
is shownbelow for the state with z =]][, 



14) ~& 



-©- 



where the solid/hollow cirlces indicate control on f / J,, and 

A f -- 



1 

e~ lhe 



(26) 



8 



The multiple-controlled-NOT gates can be constructed using 3-qubit Toffoli gates, additional ancilla (M — 1 
gates/ancilla for Mcontrol qubits) and a controlled-NOT and (see [27| page 184). 

For the hopping terms, we must simulate unitaries which implement evolution according to some product of r x 's 
and r y 's on some subset of walker qubits, if the other qubits are in some given state - a multi-qubit controlled 
operation. Firstly, the evolution by the Hamiltonian consisting of a product over r z operators can be simulated using 
controlled-NOT gates and a phase gate with a single ancilla [27j ]. 
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which outputs exp [— ifter-f t| t|] Using U exp[— iV]U^ = exp[— iUVU^] for unitaries U and V, we can use single 
qubit gates and the circuit above to simulate any product of f t x 's and r y 's. Since controlled-NOT is its own inverse, 
the controlled evolution is implemented by simply making the A e a controlled gate, i.e. 



- U 



V 
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gives exp[-ihePlF\F\T$TgT$} for Ut z W = f x and Vt z V^ = f. 

The complexity of the circuit to simulate a quantum walk will depend upon the graph, and how the nodes are 
labelled. One simplification is to minimize the Hamming weight (number of different bits) between connected nodes, 
which we use below for the walk on the line and hyperlattice. 



3. Hyperlattice walks mapped to qubits and gates 

We start with a line with 2 N nodes such that the general Hamiltonian is H = — 1 Aj[cfci+i + H.c] + tic\ci. 

The encoding of the node states is as follows: start with a single qubit, defining a two node walk, with the nodes 
labelled as ||) and This quantum walk is simply defined by H = — AiTf . Now add an additional qubit, such that 
each node now has two labels, without changing the Hamiltonian, we have two, two node walks, which we now join 
together at opposite ends, such that the order of the nodes is now J, J,, J, f, ft, and fj,. We then continue in this fashion 
(as shown in the figure ^ for N-qubits, giving a 2 N node walk on the line. Note that the label of each node differs 
from it's nearest neighbours in only one bit. Given a bit-string x = xn%n-i ■ ■ -X2X\ specifying a node, the position 

• • — ^ • • • • — ► • • • • • • • • 

II) It) lli> lit) Itt) lti> Wll) lilt) litt) liti) Itti) Ittt) Itit) Itii) 

FIG. 2: Encoding for quantum walk on the line, using 1, 2 and 3 qubits. 
along the line (with H . . . I corresponding to the origin, ie., position 1) is given by the function 

N /n-1 \ 

F(x)=l + '£2 N - n [@ XN A, (27) 

ra=l \j=0 / 

where © denotes addition modulo 2. 
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This labelling results in the following 7V-qubit Hamiltonian for the quantum walk on the line, 



E f E 

m=l \ x:F(x)=2 m - 1 (2n+l),n=0,l,... 



Af (s) t+ J] P n" +H.c. + ]Te F(s) \x)(x\ 



(28) 



such that each hopping term consists of only one Pauli term, and the rest projection operators. For the corresponding 
circuit simulation, this means that only multiply-controlled single-qubit gates arc required. In the case of uniform 
hopping, Aj = Aq, the sum over the hopping terms simplifies to 



H hop = -2 (ff + f? p} 1} + n 



x TO (2) p (D + f xp(3) p (2) p (l) + + f , p (iV-l) p ( J V-2) 



•<0 



(29) 



The corresponding circuit to simulate C4(e) = exp (— iheHk), for e = t/iV with = f^pj- fe ^pj* 2 ^ . . .P^, (such 
that the corresponding unitaries Uk are controlled rotations on the k th qubit.) is shown below (for 6 qubits); 



— X 4£ 



X4t 



x~ 



x 4e 



Xn 



X 4 e 



Xtt 



X v 



X^ — 



x w 



X 4e 



where we have used the notation Xq = R x (9) = exp(— ih0r x /2), such that X^ gate corresponds to the Pauli-X i.e. 
a bit flip. 

To write the circuit above in terms of a one- and two-qubit gates we use the construction described above. Explicitly, 
we require the multiply controlled gate 
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with the Toffoli gates realised using single qubit rotations and CNOT gates, as shown below: 



- C -9- B -0- A 



a -9- 



C -©- B- 1 A 



a -9 
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using the following single qubit gates (where R a (0) = exp(— i6a a /2)) 

A = R z (^)R y (^) 



B 




C=Ry{\) 




' 1 




' 1 


a = 


e-"/ 4 


, -a = 


e"/ 4 


R x (2e), which is 


simply: 





- x 2e - - f -e- e -e- d - 



where 



This can be simply modified to the quantum walk on the circle, by modifying the last term in the Hamiltonian to 
f^-ATP^ 2 -* . . .Pj "*, and in turn altering the corresponding gate. Having the hopping amplitudes between nodes equal 
greatly simplifies the quantum circuit simulation - the number of gates requires scales approximately as 0{n 2 ) for 
each incremental time step. 

The construction of the qubit quantum circuit for simulating the quantum walk in the line can be easily generalised 
to simulate a quantum walk on an arbitrary D-dimensional hyperlattice, with 2 nd nodes. 

Each node on the hyperlattice is specified by D bit-strings of length N, each of which denote the location of the 
node in a given direction - each node is represented by an N x D qubit state, \xi\ X2', ■ ■ ■ ; Xd), where Xk is an iV-bit 
string. 




FIG. 3: Encoding for quantum walk on the two dimensional lattice. Each node is encoded via two bit strings, of length 3 in 
this case. 



Using this encoding, the quantum walk on the hyperlattice simply corresponds to D individual quantum walks on 
the line, where D is the dimension of the lattice - there is no interaction between qubits specifying different directions. 
Thus, we use the above construction on D different sets of M-qubits to define the quantum walk on the D-dimensional 
hyperlattice as follows 



D 



H = J2 



N 

E 



E 



A 



F{x d ) T m d 



\x d :F(x d )=2" 



1 (2n+l),n=0,l, 



n 1 

n d ^m d 



H.c. 



E £F (^) \ X d)i x d\ 



(30) 



We have discussed the construction of qubit Hamiltonians for a given walk when the graph structure is completely 
known. Another scenario is where we are given access to a 'black-box' or oracle, which contains information about 
the graph structure, e.g the adjacency matrix. In the standard set-up, we may query the oracle with two nodes to 
determine if there is such a connection. This is the situation in the Childs et al. algorithm Q, and was considered 
more generally by Kendon [15[ . 
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IV. FROM QUBIT HAMILTONIANS TO QUANTUM WALKS 

The other direction to approach these mappings from is to start with a multi-qubit Hamiltonian, and determine a 
corresponding quantum walk. We begin with a simple one-dimensional spin chain. 

A. Spin-chains to Quantum Walks 

The XY model in one dimension corresponds to a chain of N qubits (spin-| particles) with nearest-neighbour 
couplings, described by the Hamiltonian 

E-J(^i+W + ^. ( 31 ) 
i=i 



which assumes homogenous coupling strengths, J. This model is exactly solvable using the Jordan- Wigner transfor- 
mation, mapping the model to a system of spinless fermions. In this representation, the Hamiltonian has a natural 
quantum walk interpretation, as fermions hopping between sites. The Jordan- Wigner transformation defines the 
fcrmionic operators 

^=(n>/)^> 4=(l[^r r , (32) 

which respect the fermionic canonical commutation relations, {fii,cj} = Sij and {£i,Cj} = {cj,cj} = 0. The spin 
operators are expressed as 

f? - I-2c\ci, (33) 

K = U.(I-^cj)ci, (34) 

j<i 

V = Y[(I-2c]cj)cl (35) 

The XY Hamiltonian then becomes 

Hxy = ^+J2- J + - M i di ( 36 ) 



j<i 



N 



describing free, spinless fermions, hopping along a 1-dimensional lattice, since the total fermion number, n — J2iLi cjc») 
is conserved. 




FIG. 4: Graph for the quantum walk given by the Hamiltonian, Hxy, with 6 sites, and three excitations. 
It is interesting to consider the same system with a higher number of excitations. In this case, the dynamics is 

N\ AT, 



restricted to a subspace with dimension D = | J = n ^ N L n y_ > where n is the number of fermions/excitations. Now 
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consider each state as encoding a node of a graph, reverting to the binary-encoding. The symmetry in the system 
results in interesting graphs for the corresponding quantum walk. For example, the N — 6,n = 3 case, where nodes 
are encoded by states of the form ITTT4I4) with three spins up, and three down, is shown in figure [4] where the two 
end states correspond to ITTTI44) an d I444TTT)- We see that this graph has a tree-like structure, leading into a cube 
in the middle. The continuous quantum walk on this graph is exactly solvable. 

• • • • • • • • • • 

J] J 2 J 3 J 4 J 5 J 4 J 3 J 2 Jj 



FIG. 5: Graph for the quantum walk given by the Hamiltonian, Hxy, with 6 sites, and three excitations, reduced to a linear 
chain. The couplings are Ji = 1, J 2 — V2, Js = 4/\/6, J a = 5/3 and J5 = 2. 

It is possible to 'collapse' such a quantum walk to a biased walk along a line [9(. This corresponds to the XF-model 
with non-homogenous coupling strengths. This is done be defining column subspaces, such that states in column 
space k, are only connected to states in column spaces k — 1 and k + 1, in terms of the corresponding graph for the 
quantum walk. Site k on the line then corresponds to an equal superposition of states in the corresponding column 
subspace. The strength of the coupling between the nodes is then determined from the Hamiltonian. Figure [5] shows 
the linear chain corresponding to to the XK-Hamiltonian with six sites, in the three excitation subspace. The two 
end nodes correspond to the states ITTT444) an d II44TTT)- 



B. Static Qubit Hamiltonians to Quantum Walks 

Now let's look at more general spin systems. A system of considerable interest, both methodological and practical, 
is the general TV-qubit Hamiltonian with time-independent couplings. As an example consider the following form: 

JV 

H = + A„r*) - £ Xtf ff + V±t?t* + Viffff. (37) 

n— 1 i,j i<j 

We have not included all possible interaction terms V^ iffj here, because the algebra then becomes rather messy, 
but instead just all the terms representing different kinds of interaction: the longitudinal and transverse diagonal 
couplings Vjj and , and a representative non-diagonal Xij ■ 

It is intuitively useful, before giving the general results, to first consider just three qubits. Using the binary 
expansion encoding, where the state \k) represents the k th node on some graph, we have 

H = [(X21 + X3i + Ai) |0)(4| + (X21 - X3i + Ai) |1)(5| + (xai - X21 + Ai) |2)(6| + (A x - X 2i - Xsi) |3)(7| 
+&12 + X32 + A 2 ) |0)(2| + (X12 - X32 + A 2 ) |1)(3| + (X32 - X12 + A 2 ) |4)(6| + (A 2 - X32 - X12) |5)(7| 

+ (X18 + X23 + A3) |0)(1| + (X13 - X23 + A3) |2)(3| + (X23 - Xl3 + A3) |4)(5| + (A3 - X23 - Xis) |6)(7| + H.C.] 

+ [VA (|0)(6| + |1)(7| + |2)(4| + |3)(5|) + V 2 \ (|0)(3| + |1)(2| + |4)(7| + |5)(6|) 
+^i(|0)(5| + |l)(4| + |2)(7| + |3)(6|)+H.c] 

+ [(ei + 62 + e 3 + 4 + V} 3 + Vl) |0)(0| + (V}[ + V} 3 + 4 - ex - e 2 - e 3 ) |7)(7| 

+ (ei + e 2 - e 3 + 4 - 4 - 4) |1)(1| + (4 - 4 - 4 - £1 - e 2 + e 3 ) |6)(6| 

(ei -e 2 +e 3 - 4 + 4 - 4) |2)(2| + (-e, + e 2 - e 3 - 4 + 4 - 4) |5)<5| 

+ (ei - e 2 - e 3 - 4 " 4 + 4) |3><3j + {-e x + e 2 + e 3 - 4 - 4 + 4) |4)(4| ]. 

which is a quantum walk over a cubic lattice, with the addition of the diagonal connections, on the faces, as well as 
on-site potentials, as shown in figure[6l If we generalise now to an iV-qubit Hamiltonian of the form above, we have a 
quantum walk on a hypercube, with the addition of next-nearest neighbour connections, where the nodes are encoded 
as described earlier for the hypercube. We can re-express the Hamiltonian in the general quantum walk form as 

2 N 

H = -J2 ^ if a + 5 4) + E e i~ c h > ( 38 ) 

(ij) 3=0 




FIG. 6: Graph for the quantum walk given by the Hamiltonian (|41l) . The nodes are labelled as in figure [T] The diagonal edges 
correspond to the two-qubit terms in the Hamiltonian, while the self-loops come from the ff terms. 



where the coefficients are defined as follows. Consider the binary representation of each of the nodes, i.e. i = i\i2 ■ ■ ■ in, 
j = Jij2---j/v where i a ,jb = 0,1 corresponding to spin-up and spin-down in the qubit representation. Then, for 
1 < a,b< N, 

{A a + Z) c (— l) J ' c Xca if i a ^ ja and ib = jb Mb^ a 

V a \ if i a ^ ja and i b ^ j b and j c =i c Mc^ a,b (39) 

otherwise 

and 

N 

ej=£(-l) J '-ea + ^- i y a+j " V H ( 4 °) 

a— 1 a,fe 

The only aspect of these expressions that is not immediately obvious is the signs. 



C. Dynamic Qubit systems mapped to quantum walks 

We now consider a universal gate set in which we allow time-dependence in all the couplings. Again we do not 
consider the most general case because the results are too messy, but instead take a special case in which the qubit 
Hamiltonian has the form 

N 

h = J2 M*)^ - A i W) - E ( 41 ) 

where we have complete control over all parameters in the Hamiltonian, which are time-dependent. This is a rather 
idealised case, but will suffice for our demonstration. If every qubit is 'connected', such that there are sufficient 
coupling terms between qubits allowing entanglement between all, then the Hamiltonian is universal for quantum 
computation. The two single qubit terms allow any single-qubit unitary to be implemented, then all that is needed is 
a two-qubit entangling operation [28[ , as provided by the XX coupling. 

From this Hamiltonian, a quantum circuit will correspond to a pulse sequence, describing applications of different 
terms in the Hamiltonian. The fundamental gate set consists firstly of arbitrary x and z rotations (on the Bloch 
sphere) for each qubit, denoted 

R x {i) = cx P H 7 f72), R z {9) = exp(-ief z /2), (42) 



which can be combined to describe any single qubit unitary operation U, via 

u = J a R z (e)Rx(i)RM, 



(43) 
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for some global phase a. As well we have the two-qubit unitaries described by 



(44) 



between qubits We will construct circuits in terms of these fundamental gates, then convert the relevant pulse 
sequence into a quantum walk. 

The canonical universal gate set consists of single-qubit unitaries and the controllcd-NOT, (cnot) operation. Using 
a method from Ref. [2^ |. we show below a circuit which is equivalent to not made up gates from our fundamental 
set; 



Xk - Zjl 

2 2 



-©- 



W 



Z-» 



X- 



For compactness of notation, we set R x (9) = Xg and R z (0 = and W = V The circuit in terms 

of the fundamental gates easily becomes a pulse sequence by interpreting the angles as times of application for 
corresponding terms in the Hamiltonian. Applying R x {l) on the second qubit corresponds to switching on A2 for a 
time T such that T = — 7/2A2. When 7 is positive, we simple replace this with the angle 7' = 7 — 2n, which gives 
an equivalent rotation. Similarly, for R z (9) on the third qubit, T = 0/2e2, and for V (x) on the third and fourth 
qubits, we switch V34 on for a time T = 

We can interpret each fundamental gate in terms of a quantum walk on graph whose nodes are arranged on the 
hypercube with the specific gate determining the edges (see figure [7]). 




i>(§> 



8)1 tit > 



tl>(§} 



9- 



/itit> 



i#>( 7 ) = exp(-z 7 ff/2), V 2 i( X ) = exp^ffff), R®(0) = exp(-i^/2 



FIG. 7: Fundamental gates as variants of a quantum walk on the hypercube. 

Imagine the 2 N nodes of a quantum walk arranged on a hypercube. An R x k ^ (7) pulse switches on connections along 
edges - figure [7^1) - in a direction given by the qubit acted upon. We then have a quantum walk on this restricted 
hypercube, for a time corresponding the angle 7. 

Similarly, a V^(x) pulse 'switches on' connections along the diagonals of faces determined by the qubits acted 
upon, resulting in a different restricted quantum walk, for a time corresponding to \ ( n g ure E[2)). 

On the other hand, a R^i\o) pulse does not connect any nodes, but rather applies a relative phase to half of the 
nodes, i.e. 



R z (6)(a |0) + b |1» = e-*> |0) + be m |1» 



(45) 



This relative phase is applied to the nodes on a 'face' of the hypercube, dependent upon the qubit acted upon (see 
figure EJ3)). A quantum computation will correspond to a series of these pulses, of varying time - the analogous 
quantum walk will be over a hypercube with time-dependent edges. As an example, we consider the quantum Fourier 
transform (QFT), the essential element of Shor's factoring algorithm. 

The QFT on an orthonormal basis |0) , |1) , . . . , \N — 1) is defined by the linear operator, 



N-l 



\J) 



e ^jk/N |jy _ ^ 



(46) 



k=0 
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b'n-l) 
|jn-l) 



-| g H r 2 I |T„_i|-|T n 

|ja> 







Jn-2 




Tr>-1 







H - T2 



H -* 



(| ) +e *2vr0. 3 „ |^ 



* (|0) + e t27r0 ^"-"" |1)) 



- (jO) + e l27r0 ' J2 "' J " |1)) 

- (|0) +e' 27r0 ^"' J " |1)) 



FIG. 8: Quantum circuit for the quantum Fourier transform. At the end are n/2 SWAP gates, reordering the qubits. 
which on an arbitrary state acts as 

N-l N-l 

x i \fi ~* Y yk i fc ) ' 

j=0 



(47) 



fe=0 



where 



N-l 



Vk 



(48) 



is the (classical) discrete Fourier transform of the amplitudes Xj . This transformation is unitary, so can implemented 
on a quantum computer. 

Following the prescription from 27], to perform the QFT on a qubit quantum computer we let N = 2™, and the basis 
|0) , . . . , \N — 1) be the computation basis for n-qubits. Each j is expressed in terms of it's binary representation, 
j = ji]2 ■ ■ -j n ~ explicitly j = j\l n ~ x + j2% n ~ 2 + ■ ■■ + j n %°- We use the notation O.jkjk+i ■ ■ ■ ji to represent the 
binary fraction jU 2 + j/c+i/4 + . . . + j//2' _fe+ . This allows us two write the action of the QFT in a useful product 
representation |27| . 



|jl ■ ■ • jn) 



2"/2 



+ e 42Taj " |1)) (|0) + e l27r0 -^-^ |1)) . . . (|0) + e i27r0 ^-J« |1)) . 



(49) 



Based on this representation, an efficient circuit, shown in figure [SJ for the QFT is constructed [27|. This circuit 
utilises the Hadamard gate, H, swap gates, and controlled- Rk gates, where 



n 



1 

e i2 */ 2 " 



(50) 



We can rewrite this circuit in terms of our fundamental gate set, to derive a corresponding pulse sequence. A 
controlled- Tfc gate is given in figure [H while the SWAP gate is shown in figure [TO] 
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FIG. 9: The controlled- gate in terms of the fundamental gate set. The pulse sequence can be read directly from the circuit. 



By combining these circuits we construct the QFT circuit in terms of our fundamental gates set. This circuit can 
be interpreted as a pulse sequences, the duration of the pulses corresponding to the angles characterising the different 
gates. 

For the above example we have assumed complete control over all parameters in the Hamiltonian, with the ability 
to switch all on or off. In physical systems, this is almost surely not the case. For example, interactions may be 
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FIG. 10: The SWAP gate as a pulse sequence using our fundamental gates. 



constant, with the single qubit terms controllable. Quantum computation is still possible in this case, though pulse 
sequences will be more complicated. An interesting problem is how circuit complexity varies as further restrictions 
are placed on possible controls. The problem of constructing efficient circuits in general is a very open and active 
area of research [30j : when decoherence is included in the operation of the gates, this becomes even more interesting 
- circuits would be designed to minimise decoherence, as opposed to complexity. Naively, one would expect less gates 
to mean shorter running time and lessening the effects of decoherence. A detailed study may demonstrate that this 
is not the case. 



V. CONCLUDING REMARKS 



In this paper we have formulated quantum walks in a Hamiltonian framework, and explored the mappings that 
exist between various quantum walk systems and systems of gates and qubits. The Hamiltonian formulation possesses 
considerable advantages. We have seen that it allows a unified treatment of continuous time and discrete time 
walks, for both simple and composite quantum walk systems. It is also necessary if one wishes to make the link 
to experimental systems. This latter point becomes particularly clear when one tries to understand decoherence for 
quantum walkers, for which it is essential to set up a Hamiltonian or a Lagrangian description. 

In the paper we have concentrated on walks on hypercubes and hyperlattices. Walks on hypercubes are naturally 
mapped to systems of gates or qubits, and we have explored mappings in either direction. Walks on hyperlattices, on 
the other hand, can be mapped to qubit or gate systems, but the mappings are not so obvious - we have exhibited 
them, and thereby shown how one could construct an experimental c?-dimensional hyperlattice from a gate system. 
In the case of both hypercubes and hyperlattices we have exhibited the general methods for finding these mappings 
and their inverses, in sufficient detail that it should now be clear how to make such mappings for quantum walks on 
more general graphs. 

The practical use of our methods and results does not become completely clear until we incorporate the environment 
into our Hamiltonian description. As indicated in the introduction, this can be done in a fairly comprehensive way, by 
using a general description of environments in terms of oscillator and or spin baths. The rather lengthy results once 
this is done appear in a companion paper to this on.e[24|. Once this is done it becomes possible to solve rigourously 
for the dynamics of quantum walk systems, without using ad hoc models with external noise sources. The results can 
be pretty surprising, as shown by the results in ref.[25j for one particular example. 

Ultimately the main reason for the work in the present paper is that one can bring the work on quantum walks into 
contact with experiment, and design experimental systems able to realise different kinds of quantum walk. In parallel 
work we have done this for both a particular ion trap system, and for a particular architecture of spin qubits 31]. 
Only in this way will it be possible to fully realise the potential offered by quantum walk theory in the lab (and to 
test it experimentally!). 
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